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Abstract 

We develop the quantum inverse scattering method for the one-dimensional Hubbard 
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states are Yangian singlet. 
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1 Introduction 



The past two decades have seen a rapid development of algebraic methods for the exact 
solution of one-dimensional quantum systems. This development was partly initiated and 
most strongly influenced by the contribution of the Leningrad school, which clarified the 
fundamental meaning of the Yang-Baxter equation for the understanding of exactly solvable 
one-dimensional systems [|], ||, [3|, It culminated in the invention of quantum groups || 
which are by now generally accepted as the mathematical framework of the theory. 

The more traditional (coordinate) Bethe ansatz approach to one-dimensional quantum 
systems consists of a direct construction of eigenfunctions and yields a system of Bethe ansatz 
equations for a set of parameters which characterize these eigenfunctions and the spectrum 
of the Hamiltonian. One may ask, if the quantum inverse scattering method is more than 
an alternative way to derive the Bethe ansatz equations. In fact, as long as we are only 
interested in quantities, which are entirely determined by the spectrum of the Hamiltonian, 
we do not need the quantum inverse scattering method. The most successful attempts on 
the calculation of correlation functions (6|, [7|, however, rely heavily on the use of algebraic 
methods. 

Unfortunately, some of the physically more interesting models, like the Hubbard model ||, 
the Kondo model or the Anderson model || have only partly been capable by an 

algebraic treatment so far. In the present article we will report on some recent progress 
concerning the Hubbard model. 

The Hubbard model was solved by (coordinate) Bethe ansatz by Lieb and Wu [|ll], |i~2^, [llf|. 
The basic tools of the quantum inverse scattering method for the system under periodic 



boundary conditions, i?-matrix and monodromy matrix, were constructed by Shastry [14, 



15, [L(|] and by Olmedilla et al. [jD], 15, [L9||. Since there exists a so-called pseudo vacuum 
state, on which the monodromy matrix is acting tridiagonally and which is an eigenstate of 
its diagonal elements, an algebraic construction of eigenstates, usually called algebraic Bethe 
ansatz, should be possible according to common belief. However, algebraic Bethe ansatz, if 
we understand it in this broad sense, is not really a method. It merely means to use commu- 
tation relations between the matrix elements of the monodromy matrix, whose interpretation 
in physical terms is moreover not a priori clear, to construct the eigenstates of some suitably 
chosen generating function of a family of commuting operators, which contains the Hamilto- 



3 



nian of the system. The most sophisticated variant of algebraic Bethe ansatz was developed 
by Tarasov p0[ , who diagonalized the transfer matrix of the Izergin-Korepin model |21]. This 
work was recently generalized to the case of the Hubbard model in a remarkable article of 



Ramos and Martins [22]. Besides the fact that their proof of the "cancellation of unwanted 
terms" seems to be incomplete, we still feel unsatisfied about two points. First, the expres- 
sions for the eigenstates are of complicated recursive nature. It seems to be unlikely that 
they can be used in the calculation of correlation functions. Second, due to the complicated 
way in which the various operators contained in the monodromy matrix enter the expression 
for the eigenstates, an intuitive physical interpretation is difficult. 



Therefore we follow a different route here |23j which is based on the work of one of us on 
the fermionic nonlinear Schrodinger model ]^4|, ^jJ. We take the thermodynamic limit first 
and construct the eigenstates afterwards. This is the original idea of the quantum inverse 
scattering method |26|, ^7], [28| , which was designed in close analogy to the inverse scattering 
theory for the solution of classical integrable systems |29| . The disadvantage of this method 
is that it is so far restricted to uncorrelated vacua (ground states). It is therefore neither 
capable of relativistic models nor of non-relativistic models with general nonzero density of 
particles in the ground state. On the other hand, there are lots of advantages which make 
it highly desirable to generalize the method. We get rid of the complicated Bethe ansatz 
equations. The algebra of the elements of the monodromy matrix simplifies, and we obtain 
an intuitive interpretation of these operators. Furthermore, there is a realistic perspective to 
calculate correlation functions by use of the quantum Gelfand-Levitan equation [pOf . 



In the following section we summarize several results for the finite periodic system [17, 31], 
which will be needed later. In section 3 we describe the passage to the infinite interval. Some 



of the technical details are shifted to [Appendix A . We obtain a renormalized monodromy ma- 
trix and the commutation relations among its entries, which are encoded in a new simplified 
i?-matrix. In section 4 we identify a suitable generating function of commuting operators, 
namely the quantum determinant of a submatrix A(\) of the monodromy matrix. The com- 
mutation relations between the entries of this submatrix decouple from the rest of the algebra 
and can be written in form of an exchange relation which is generated by a submatrix r(A,/x) 
of the new i?-matrix. r(A, fi) is 4 x 4 and after an appropriate reparametrization turns into 
the rational ii-matrix of the isotropic Heisenberg spin chain. Thus ^4(A) provides a represen- 
tation of the Yangian quantum group Y(su(2)). In section 5 we identify some of the elements 
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of the monodromy matrix as creation and annihilation operators of eigenstates of the Hub- 
bard Hamiltonian. We discuss the structure of higher conserved quantities, and we calculate 
the action of the Yangian on the eigenstates. In the first part of section 6 we propose two 
pairs of normalized creation and annihilation operators of scattering states. These operators 
constitute representations of the right and left Zamolodchikov-Faddeev algebra, respectively, 
and mutually anticommute. We interpret them as generators of fermionic quasi-particles. 
The Zamolodchikov-Faddeev algebra provides their S'-matrix. These operators create the 
known Bethe ansatz states on the infinite interval. We work out the action of the Yangian 
on the Bethe ansatz states. The Yangian mixes spin multiplets, which are degenerate in the 
thermodynamic limit. In the second part of section 6 we propose creation and annihilation 
operators of bound states of quasi-particles. These bound states correspond to the string 
states of coordinate Bethe ansatz |^2|. We calculate their S'-matrix. All of them are Yan- 
gian singlet. Section 7 is left for a summary and outlook. Appendix B contains a list of 
all commutation relations between the elements of the monodromy matrix. Appendix C is 
devoted to a discussion of the construction of bound state operators. In Appendix D we give 
explicit expressions for some higher conserved operators, which are needed for the discussion 
in section 5. 



2 Hamiltonian and Monodromy Matrix under Periodic Bound- 
ary Conditions 

The Hubbard model in its most basic single-band version is describing the dynamics of inter- 
acting electrons inside the conduction band of a solid. Its Hamiltonian is usually formulated 
in terms of creation and annihilation operators Cj a , Cj a of electrons in Wannier states, 



j>=T4 o 



1\ / 1\ 1 



n n-- 2 \nn- 2j 4 



(2.1) 



The index j runs over all Wannier states which may be identified with the lattice sites of the 
solid, a =|, J, is the spin index and rij a = c^ a Cj a is the density operator. The interaction part 
of the Hamiltonian is thought to model the screened Coulomb interaction of the electrons. In 
the following we will use the terminology of quantum field theory, and we will call the empty 
band state |0) the zero density vacuum. Since we want to study finite excitations over the 
zero density vacuum |0) in the thermodynamic limit, we normalized the Hamiltonian such 
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that H\Q) = 0. 



The most distinctive feature of the one-dimensional model is that its Hamiltonian may be 
embedded into a family of mutually commuting operators, which is generated by a properly 
constructed transfer matrix. This feature is commonly referred to as quantum integrability. 
Below we will give a short account of the construction of the transfer matrix and its related 
monodromy matrix for the Hubbard model under periodic boundary conditions. We start 
from the local exchange relation |l£ 

ft(A,M)[£j(A)® s £i(M)] = [£ i (//)® s £ i (A)]ft(A,M). (2.2) 
The symbol <8> s in this equation denotes the Grassmann direct product 

[A ®s B] aj>f3S = (-l)^ +p ^ p ^A aP B l5 (2.3) 



with grading P(l) = P(4) = 0, P(2) = P(3) = 1. The P-matrix 1Z is a c-number matrix 
which encodes the commutation rules of the matrix elements of the L-matrix Cj. These 
matrix elements are operators acting on the j-th Wannier state. We adopt the expressions 
for the matrices 1Z and Cj in terms of two parameterizing functions a(A), 7(A) from ref. 18]. 
For later convenience, however, we shift the arguments of a(X) and 7(A) by ?, such that we 
simply have a (A) = cos A, 7(A) = sin A. Then the L-matrix is 



where jjv(A) = (1 
defined as 



_ie^(A) c t J 



e " :X: fr<i, 

p -h(X)J 



t 



lc iT c il e 

-9nm 



gh(A) 



(2.4) 



j(j) sin A + irijo- cos A, <7j CT (A) = (1 — n JO -) cos A — in^ sin A, and h(X) is 



sinh2/i(A) _ U 
sin 2A ~ T ' 



(2.5) 



Due to space limitations we do not reproduce the P-matrix here. It is 16 x 16 and contains 
36 non-vanishing entries, only ten of which are different modulo signs. The ten different 
entries are denoted by pi, i = 1, ■ • • , 10, in ref. |18[| . They are rational functions of cos A, sin A 
and e h ( x \ A list of the matrix elements and some basic formulae which have been used in 



our calculations can be found in Appendix A of ref. [23]. Eq.( |2.2| ) considered as an abstract 
definition of an algebra has the property that tensor products of representation of this algebra 
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are again representations. This property is called co-multiplication property. It assures that 
the monodromy matrix 

r mn (A) = £ m _i(A)£ m - 2 (A)---£ n (A) (ro > n) (2.6) 

satisfies again eq.(|2.2j). T mn {\) contains all information about the Hubbard model under pe- 
riodic boundary conditions. After proper renormalization T mn {\) turns into the monodromy 
matrix of the Hubbard model in the thermodynamic limit, which will be the central ob- 
ject of investigation of the present paper. Before we continue with the description of the 



thermodynamic limit, we list the most important properties of T mn (X) [31]. 
The transfer matrix 

r mn (A) = str(T mn (A)) = tr((a z <g> a z )T mn {X)) (2.7) 

generates a family of mutually commuting operators M, [l^, [l8|, [3l]] , 

if a 

mr mn (A) = —(N mn -2(m-n)) + iU mn 

+ \(H mn + {m-n) l £j +0(A 2 ). (2.8) 

Here 

m— 1 

N mn = ( n ;t + n n) (2-9) 

j=n 

is the particle number operator, and Tl mn is the lattice momentum operator. For the subtle 
question how to define II mn properly such that it commutes with the Hamiltonian we refer 



the reader to Ref. [31]. H mn is the Hamiltonian (|2.1[) on a one-dimensional lattice of length 



m — n under periodic boundary conditions (let j run from n to m — 1 in eq. (|2.1| ) and let 



Clearly, the Hamiltonian (^^) is invariant under su(2)-rotations generated by the opera- 
tors of total spin 

1 m— 1 

sa = 2 E <A^- ( 2 - 10 ) 

j=n 

Here a = x,y,z and the matrices a a are the Pauli matrices. Due to the invariance of the 
Hamiltonian under the transformation 



■#t-^ T> c n ^{-\yc\ v U^-U (2.11) 
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there exists a second su(2) symmetry |33|, [34|, |35|. This symmetry is called repairing symmetry. 
Applying Q2.ll ) to ( |2.10[ ) we get its generators if in the form 

-, m— 1 

n* = .E(-i) i (^| l + ( i-i c it). ( 2 - 12 ) 



2 



'/" - lEt-^W-^t)- (2-13) 



■ m— 1 

2 

r/ 2 = ^(N mn -m + n). (2.14) 



Note that the transformation ( p. 11 ) twists the boundary conditions, ifm — n is odd. Therefore 



H mn commutes with rf and rf only if the lattice has an even number of sites. Since rf is 
essentially the particle number operator, we may understand the ry-pairing symmetry as a 
non-Abelian extension of gauge symmetry. 

Rotational symmetry and ry-pairing symmetry both extend to symmetries of the mon- 



odromy matrix [31]. In order to make this statement explicit, we have to introduce certain 
matrix representations of su(2). We will denote the n x n unit matrix by I n . Let 

E x s = ^{a + ®a- + <j- ®a + ), (2.15) 
= - l -{a + ®a- -a' ®<t+), (2.16) 
1 (a z ®I 2 -I 2 ®(T Z ), (2.17) 



and 
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= ~(<7+ ®a + +<j- (2.18) 



2 

2 

4 

These matrices obviously satisfy the su(2) commutation rules 



= -Ua + ® <t + - a~ ® a~), (2.19) 
= \{a z ®/ 2 + / 2 ®<7*)- ( 2 - 2 °) 



[£?,£?] =iE ate E? J j = S ,r/. (2.21) 



Let us perform a basis transformation 



£a = G afc s fe ^ 222 ^ 
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with transformation matrix 



G 



0-10 
1 
1 



(2.23) 



Since G is orthogonal with det(G) = 1, the transformed matrices S" satisfy (2.21). We are 
now able to state the su(2)©su(2) symmetry of the monodromy matrix T mn (X), 



[T mn (X),t a s +S a }=0, 
[T mn (\),t a v + v a }=0. 



(2.24) 
(2.25) 



If a = x,y in ( 2.25] ) we have to require both m and n to be odd. ( |2.24| ) and (2.25) imply 
the invariance of all higher conserved quantities in the expansion ( [2.8D under rotations and 
repairing transformations. 



The twist ( ggg ) may appear somewhat unnatural. However, we had to introduce it here, 
since we wanted to keep the notation of the earlier paper [23|. We may remove the twist by 
a gauge transformation in auxiliary space. Let 



W = diag(l,l,i,i). 



Then 



(2.26) 



(2.27) 



Thus the gauge transformed monodromy matrix WT mn (X)W 1 satisfies ( p. 24 ) and fl2,25|) 



with £" replaced by Note that W <S> W commutes with the .R-matrix |31j| , which implies 
that the exchange relations for T mn (X) and WT mn (X)W~ 1 are the same. 

The grading of the monodromy matrix, its behavior under particle-hole transforma- 
tions [31] and the structure of the matrices suggest the following block notation for the 
monodromy matrix, 

( D n (X) C n (A) C 12 (A) L» 12 (A) \, 

flu (A) An (A) A 12 (X) B 12 (X) 

B 2l {X) A 2l (X) A 22 (X) B 22 {X) 

\ D 21 (X) C 2 i(A) C 22 (A) D 22 {X) ) 

We will see in the following that many of the algebraic properties of the Hubbard model are 
conveniently expressed in terms of the 2x2 submatrices -A(A), B(X), C(A), -D(A). As an 



^mn(A) 



(2.28) 
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example let us describe the behavior of T mn (X) under hermitian conjugation, which will be 
needed later and which can be obtained by using the methods outlined in Ref. |31], 



(5(A))t 
(C(A))t 
P(A)) t 



a y A(w/2 - X*)a v , 
ia y B{Tr/2 - X*)a y , 
ia y C(vr/2 - X*)a y , 
a y D(n/2 - X*)a v . 



(2.29) 
(2.30) 
(2.31) 
(2.32) 



The dagger on the lhs of these equations means hermitian conjugation in quantum space, 
and the asterisk on the rhs denotes complex conjugation. For notational convenience we did 



not attach labels m and n to the submatrices A(X), ■ ■ ■ ,D(X) on the rhs of (2.28). We will 
keep the same notation below, when we discuss the thermodynamic limit. 



3 Passage to the Infinite Interval 

As we shall see in the sequel, carrying out the thermodynamic limit leads to a severe simpli- 
fication of the i?-matrix. The commutation relations between the entries of the monodromy 
matrix will become simple enough to allow for an identification of creation and annihila- 
tion operators of quasi-particles, generators of conserved quantities and symmetry operators. 
The thermodynamic limit cannot be taken naively. The monodromy matrix requires infrared 
renormalization, which has to be done with respect to a given vacuum characterized by macro- 
scopic parameters. These are the density of electrons pn and the magnetization density pm- 
As result of the thermodynamic limit we will obtain the finite energy excitations over the 
chosen vacuum. In contrast to the algebraic Bethe ansatz for the finite periodic system we 
will not be able anymore to distinguish between a pseudo- vacuum, upon which all eigenstates 
of the transfer matrix are built by the action of creation operators, and the physical vacuum, 
which is the true ground state of the model. In general, in the thermodynamic limit both 
states will be characterized by different values of pm, Pn and thus will be separated by an 
infinite energy difference. 

Infrared renormalization of the monodromy matrix is done in analogy with the inverse 
scattering method for integrable classical systems |29| by splitting off the asymptotics of the 
vacuum expectation value of the monodromy matrix for m, —n — > oo, which therefore has 



to be known a priori. For this reason the method [26, g7], £g] is so far restricted to 
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(asymptotically) uncorrelated vacua. In case of the Hubbard model there are four possible 
choices, the empty band (pu = Pn = 0), the completely filled band (pm = 0,Pn = 2), and 
the half-filled band with all spins up (pm = 1> PN = 1) or an spins down (pm = —1, Pn = !)■ 
In the following we will restrict ourselves to the empty band vacuum |0) which is defined by 

c ja \Q)= 0. (3.1) 



Our description of the general method closely follows Sklyanin |28| . Define the Hilbert 
space 7i of states of "compact support" as the space of all finite linear combinations of vectors 
c mi<ri ' ' ' c ln n a n Denote the vacuum expectation value of the L-matrix by 

V{X) = <0|£ m (A)|0). (3.2) 

V(X) does not depend on m, because of the translational invariance of the vacuum. 
Let 

4(A) = vixy^q^vixy, (3.3) 

Znn(\) = V(\r m T mn (X)V(X) n . (3.4) 

It is easy to see that the limits lim n ^_ 00 (j;|7^ tri (A)|?/) and lim rrt ^ 00 (x|7^ in (A)|?/) exist for all 
\x), \y) € TL. These weak limits determine a pair of operators 

t+(A) = lim T nm (A), (3.5) 

n— >+oo 

f~{X) = lim f mn (X), (3.6) 

n— )— 00 

with asymptotics 

lim t+(A)= lim f m (X) = h. (3.7) 

Multiplying ( |2.6| ) from the left by C m (X) or from the right by £ n _i(A), respectively, we obtain 
two recursion relations for T mn (X), which imply a pair of recursion relations for (A) and 



T m (A). By use of the asymptotic condition (|3.7|) these are equivalent to the following pair of 
Volterra "integral equations" for 7^ (A), 

oo 

t+(A) = / 4 + E ^ + ( A ) (A-i(A)-I 4 ), (3.8) 

j'=m+l 
m— 1 

T-(A) = J 4 + E (4(A) -A) V(A). (3.9) 

j=-oo 
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The above considerations imply the existence of the weak limit 

T(A) = lim T mn (A) = r+(A)T-(A). 



m,— n-^roo 



(3.10) 



T(A) is the renormalized monodromy matrix. Equation ( |3.8[) , or (|3.9| ) respectively, implies 
the "integral representation" 



f (A) = / 4 + £(£ m (A) - I 4 ) + E (^™( A ) - h)(C n (X) - h) + ■ ■ ■ ■ 

m m>n 

Note that (0|(£ m (A) - J 4 )|0) = by construction. Hence (0|T(A)|0) = h- 
£ m (A) can be easily calculated. We find 

V(X) = diag(-e h(A) sin 2 A, e _/l(A) cos A sin A, e h ^ cos A sin A,-e h(A) cos 2 A), 
and thus 

An(A) = F(A)~ m - 1 /: m (A)y(A) m 



(i cot A) 



,t e^W -imp(A) 



-i(icotAr-T C ^|^ e 
4 T (icotAr^Se- imp(A) 



cos A 

V ^T^id-A^A 



(icotA^Tc^^e^M 
(icot A) nm T^ n ™.i 

^-m^rnl s i n AcosA 



-ic mT (icotA)"-i^ r e 

^THnJ. sinAcosA 



4 T (icotA)-^i^e- imfe W 

n m j e~ fa ( A ) imp(A) _■ 1 rn< -2m \ \ 



ic mT (icot A)-^f^e imfc W 
(icotA)-^Tc mi |^e imfc ( A ) 



-i(icotA)- n -T c J n| ^ A le- imfc ( A ) (icot A)-""*""™* ^ 



(3.11) 



(3.12) 



(3.13) 



(icot A)- n ™T +n ™J- 

ml cos A 

Here we have introduced new functions 

e ife(A) = _ e 2h(X) CQt Aj e i P (A) = _ e -2fc(A) CQt Aj (3 14 ) 

which we adopted from the recent analytic Bethe Ansatz for the Hubbard model by Yue and 
Deguchi g§. 

Now we will turn to the calculation of the commutation relations between the elements 
off (A). Let 

Cg\\,n) = £ m (A)® s £ m Gu), (3.15) 
(A, a*) = T mn (A)® s T mn (/x). (3.16) 
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(2) 

We may apply the above discussion of the renormalization of T mn {X) to Tmn(X,fi), if we 
replace V(X) by 

V^^fi) = (0|£ W (A) ® s £ m (/i)|0}. (3.17) 

Note that V^ 2 \X,fi) is not just the tensor product V(A) (x) s V(//). There appear additional 
off-diagonal terms due to normal ordering of the operators. We obtain a renormalized tensor 
product matrix 

fW{\^)= lim V {2 \\^)- m T^{\^)V^{X^) n , (3.18) 

m,—n—>oo 

which satisfies (0|T( 2 )(A, /i)|0) = h 6 . Taking the vacuum expectation value of the local 
exchange relation (|2.2D yields 



n{\^)v^{\^) = v^i^xmx,^, (3.19) 

and we conclude that 

K{\,(j,)f( 2 \\,») = f®(v,,\)K{\,(J,). (3.20) 

If 7mn(A, //) is defined in analogy to T mn (X) with V(A) replaced by V^ 2 \\, fx) in the definition 
flU), then 

t mn (A) ® s T mn (/x) = C/ m (A, A i)- 1 t£)(A,^)?7 ri (A, At ), (3.21) 

where we have set 

£/ n (A, M ) = V^(X^)- n [V(X) n ® s (3.22) 
Assume for a while that the limits 

U + (X,/j,)~ 1 = lim U m (X, /j,)" 1 , f/-(A,/x)= lim U m (X,n). (3.23) 



m— >— oo 



exist in a common domain of convergence. Then, according to eq. ( |3.21| ) T mn (X) ® s T mn (n) 
has a weak limit for m, —n — ► oo. We identify this limit with T(A) (8> s T(n), 

f(X) ® b T(m) = t/+(A, / i)- 1 T( 2 )(A, At )?7_(A, At ). (3.24) 



Finally, inserting the above equation into ( p. 20 ), we arrive at the exchange relation for the 



monodromy matrix 7~(A) on the infinite interval, 

£ (+) (A,/i) fr(A)® 8 r(M)l = [r(/x)® s r(A)l^-)(A,/x), (3.25) 



where 



ftM(A, /i) = C/±(/x, A)-^(A, fi)U±(\, n). 



(3.26) 
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The calculation of the matrices U±(X, fi) is rather technical. We present it in Appendix A 



Here we only note three important facts, (i) there is no common domain of convergence for 
all matrix elements of C/+(A, fi) and U-(X, //). We will come back to this point later, (ii) if we 
stay away from some singular points (cf. [Appendix A ) then C/+(A, p) a p^8 = U-(X, fJ,) a /3,'y6- 
(hi) it is a nontrivial matter of fact that all the matrix elements of U±(X, p) are simple rational 
functions of the original Boltzmann weights pj{X,jj). 



Using the explicit form of the matrices U±(X, p,) given in Appendix A , we obtain 

n(\n) = n {+ \x,p) = tz^-\x^) = 



I PI 












































o\ 














P1P4 
ipio 
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v° 












































pj 



(3.27) 



The reader is urged to compare this expression with the i?-matrix on the finite interval 18]. 
Instead of the 36 non-vanishing elements of the original i?-matrix we have only 18 non- 
vanishing elements here, which brings about simpler commutation relations between the 
elements of the monodromy matrix. All matrix elements except the two diagonal elements 

2 

p ' ip4 _ p2 are just at the position of the l's of the permutation matrix, which means that, 

2 

were it not for the two elements p3pA p2 , all commutation relations would reduce to the 

P3-P1 ' 

mere interchange of two factors along with a multiplication by some rational function of the 
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Boltzmann weights. A close investigation of (|3.25| ) shows in particular, that the elements 
of the submatrices A(X), ■ ■ ■ , D(X) generate sub-algebras of the exchange relation. The sub- 
algebras generated by A(X) and by -D(A) are again Yang-Baxter algebras with certain new 
i?-matrices which are submatrices of TZ{X,(j). The sub-algebras generated by B(X) and 
C(A) can, after certain normalizations, be identified as representations of the right and left 
Zamolodchikov-Faddeev algebra. The construction of these algebras and the discussion of 
their physical meaning will be the contents of the following sections. 

Let us come back to the remark (i) above. Due to the peculiar convergence properties of 
the limits U±(\, /j.) the exchange relation in the form ( 3.25[) has only symbolical meaning. It 



has to be understood as generating a set of 256 equations, each of which is meaningful. These 
equations, however, are not necessarily defined on the same domain in the A, [i parameter 
space (cf. |Appendix A| ). 

4 Yangian Symmetry and Commuting Operators 

The definite goal of this work is to construct algebraically the eigenstates of the Hubbard 
Hamiltonian (|2.1|). As usual in the theory of integrable systems, we will not directly work 
with the Hamiltonian, but with an appropriately chosen generating function of a whole family 
of mutually commuting operators. For the finite periodic system this generating function is 
the logarithm of the transfer matrix r mn (A). Commuting operators are obtained as the 
coefficients of its expansion around A = fl2.8|) . Equations ( |3.11|) and ( |3.13| ) show that this 
expansion does not exist for the renormalized monodromy matrix T(X). There is substitute, 
however, which is intimately connected with the existence of an additional Y(su(2)) quantum 



group symmetry of the Hubbard model in the thermodynamic limit [23]. We shall describe 
it below. 

The commutation relations between the elements of the submatrix A(X) decouple from 
the rest of the algebra. 

r(A, M ) (A(X) ® A(//)) = (A(fji) ® A (A)) r(A, fi), (4.1) 



where 



(4.2) 
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and V is a 4 x 4 permutation matrix (Vx <g> y = y®x). If we introduce the reparametrization 



V (A) = -2i cot 2A cosh 2/i(A) 

= -2 sin /c(A) + y = -2 sinp(A) - y , 



(4.3) 



the ii-matrix r(A, //) turns into the rational i?-matrix of the isotropic Heisenberg spin chain, 

i(7 + v (A) - u( / u) 

Let us assume that ^4(A) allows for the following asymptotic expansion in terms of v (A), 



A(X) =I 2 + iUj2 



1 



n+l 



(4.5) 



n=0 ^ ' \a=l / 

where a a = G ab a b (cf. (^23|)). Then it follows from general considerations p7| , pB] , Q| that 
the first six operators Qo,Qf generate a representation of the Y(su(2)) Yangian quantum 
group. 



There is the following alternative description of the Yangian Y(su(2)) f39|| . The Yangian 
Y(su(2)) is a Hopf algebra which is spanned by six generators Q° (n = 0, 1, a = x,y,z), 
satisfying the following relations, 



V0' HQ — J V0' 

Q a o,Q\] =f abc Q c i, 
QiiQ\ , Qo,Q\ 



+ 



Qi,Qi , Qq,Qi 



K 2^ A abkefgjcdk + jcdkefg jabk} |ge ^ Q^Qdy 



(4.6) 

(4.7) 
(4.8) 



Here n is a nonzero constant, f abc = ie abc is the antisymmetric tensor of structure constants 
of su(2), and A abcde f = f adk f hel f°f m f klm , The bracket { } in (fl|) denotes the symmetrized 
product 

! ^ (4-9) 



{XI,X 2 ,X 3 } = — ^ X a lX a2 Xa3- 



Being a Hopf algebra Y(su(2)) carries an outer structure (co-multiplication, antipode, co- 
unit), which is described in ref. |39| and which assures that Y(su(2)) has a rich representation 
theory |5| gl|. 

A careful consideration of the limit v(\) — > oo shows that ^4(A) is indeed of asymptotic 
form ( |4.5| ). There are several possibilities to carry out this limit as a function of A. We 
found, however, that only one of these yields finite results for Qq and Qf. We have to take 
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9(A) — > oo and have to choose the proper branch of solution of eq. ( |2.5|) , which determines 
h as a function of A. (2.5) implies that 



-2h(X) 



— sin2A ±Jl + (jsin2A 



(4.10) 



in order to achieve convergence of the matrix elements T a p (a, j3 = 2, 3) we have to choose 
the lower sign here, then e~ 2h W is approximately equal to — ^ sin2A for large positive values 



of sin2A, and we obtain 



a 2h(A) 



| e 2iA + 0(e 6iA ); e ifc(A) = _i {e 2iA + 2e 4iA + 0(e 6iA )}; 



e -m = ±{e 2 ^ - 2e 4iA + 0(e 6iA )}, ^ 



-^{e 2iA + 0(e 6iA )}. 



(4.11) 



The leading terms in the series (3.11) are of order e 2lA , e 4lA , • • •. Thus, from the first two sums 
in ( |3.11 ), we get the expansion of the matrix A(A) up to order e 4lA , and the last equation in 
(|4.11 ) yields the required expansion in (v(X))^ 1 up to second order. 



The final result for the representation of Yangian generators is 



Qi 



l 



(4.12) 



3 i>3 

The factor \U occurring in ( [4.13 ) can be identified with the constant k in ( |4.8| ). Note that 
Qq = S a is just the operator of the a-component of the total spin (cf. ( |2.10| )). The Yangian 
representation ( [4.12 ) and ( 4.13 ) was first obtained by Uglov and Korepin |)42| , |23|| . It can be 
embedded into a larger family of Yangian representations connected with long-range-hopping 
extensions of the Hamiltonian ( |2 . 1|) |43|| . Uglov and Korepin showed that Qq and Qf commute 
with the Hamiltonian on the infinite line. 



Since the quantum determinant 

DeVKA) = A n (X)A 22 (X) - A 12 (X)A 21 (X), 
where v(X) = v(X) — ill, is in the center of the Yangian 

[Det q A(X),A(fj,)] = 0, 
and thus provides a generating function of mutually commuting operators, 

[Det q A(X),Det q A(n)] = 0, 



(4.14) 



(4.15) 



(4.16) 
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it is a natural candidate to generate the Hamiltonian, too. Performing again the asymptotic 
expansion in terms of v(\)~ 1 , 

CO j 

DeV4(A) = l + ify£— ^ TT , (4.17) 

we obtain Jo = 0, J\ = iH, i.e. the Hamiltonian is indeed among the commuting operators 
generated by T)et q A(X) . All the conserved operators are Yangian invariant by construction. 
We discuss their relation to the formerly known conserved quantities 1 1£ , [ll| 44, 45 1 in section 
5.2 below. 

In closing this section we shall add a comment. The Hubbard Hamiltonian on the infi- 
nite interval is invariant (up to a constant) under the transformation ( 2.11|) . The Yangian 



generators Q$ and Qf, however, are transformed into a pair of generators Qq and Qf of 
a second, independent representation of Y(su(2)) p2|| . These two representations mutually 
commute. Therefore they can be combined to a direct sum Y(su(2))©Y(su(2)). The reason 
why we get only one of these representations from our QISM approach is that, in order to 
perform the passage to the infinite interval, we refer to the zero density vacuum |0). This 
vacuum has lower symmetry than the Hamiltonian. It is invariant under the su(2) Lie algebra 
of rotations, but does not respect the ry-pairing su(2) symmetry of the Hamiltonian. A fully 



su(2)©su(2) invariant vacuum would be the singlet ground state at half filling [46]. It seems 
to be yet a formidable task to formulate the QISM with respect to this state. 



5 Conserved Quantities and Eigenvectors 
5.1 Eigenvectors of the Hamiltonian 

In the present section we will try to understand the meaning of the operators contained in 
B(X),C(X), D(X). To begin with, let us have a look at the commutator of T(X) with the 
particle number operator N, which is the extension of N nm (|2.9| ) to the infinite interval. 
Note that [V(X), E*] = 0, and thus by ( ^25|) 

[N,f(X))=2[f(X),m (5.1) 

This is a set of 16 equations. Writing out the equations explicitly, we find that -Dn(A), 
-022(A) and A{X) conserve the number of particles. The operators B a \{X) and C2a{X) increase 
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the number of particles by one, whereas B a 2(X) and Ci (A) reduce the number of particles 
by one. L>2i(A) adds two particles to the system, whereas Di2(A) removes two particles. 
Hence, B a i(X),C2 a (X) and -D2i(A) are candidates for creation operators of eigenstates of the 
Hamiltonian. The action of these operators on the vacuum follows from ( |3.11| ), 
ie fc(A) 



Bll(A)|0) = -^sT^ e " imp(A)c ^ |0) ' (5 ' 2) 

li(A) 

^2i(A)|0) = ^J2 e ~ imPiX)c U I0>, (5.3) 



1)1 

MX) 

ooTA^ ^ 

C 21 (X)\0) = - e ^~ imm 4n^ (5-4) 

m 

C 22 (A)|0) = ^"^W^JO), (5.5) 

£>2i(A)|0> = ^2^ E c lAl { 6{ - m ^ n)e-' l( - mk M +np W + 6(m < n)e~ i( - nk ^ +mp ^j^) 

These states are the most elementary eigenstates of the Hamiltonian as can be verified by 
direct calculation. Note that they are bounded in different, disconnected parts of the spectral 
parameter space. This is because of the constraint (|2.5|), which turns into 

sin/c(A) — sinp(A) = X — (5-7) 

when expressed in terms of k(X), p(X) ( |3,14| ). If k(X) is real, p(X) cannot be real for non-zero 
U and vice versa. D2i(X) creates a bound state, if Q(k) = —$s(p) < 0. This condition is 
compatible with the constraint ( |5.7|) on k(X) and p(X). Let U ^ 0. Then k := $s(p) = —$s(k) 
implies q := 3?(/c) = 3?(p)mod2-7r, and 

sinliK = — . (5.8) 

4 cos q 

Hence, if U < 0, D2i(X) creates a bound state for \q\ < tt/2, and if U > 0, D2i(X) creates a 
bound state for tt/2 < \q\ < tt. 

The commutators of the various operators contained in 7~(A) with Det g (A(/x)) are sum- 
marized in Appendix B. Let us insert the asymptotic expansion ( [4.17 ) into ( B.23| )-( jB.28 ). 
Then there occur only two different ratios of Boltzmann weights in these equations, which 
may be expanded by use of ( |4.11| ) as 

-m - 1+ ^y+(f -*-™)^°(^ ^ 
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Comparing terms of second order in v(fi) , we obtain the following commutators, 

[H,B al X)] = -(2cosp(A) + *7/2)£ ol (A), (5.11) 

[H, B a2 (X)} = (2cosk(X) + U/2)B a2 (X), (5.12) 

[H,C la X)} = (2cosp(X) + U/2)C la (X), (5.13) 

[H,C 2a (X)} = -(2cosk(X) + U/2)C 2a (X), (5.14) 

[H,D 12 (X)} = 2(e^+e~ ik ^)D 12 (X), (5.15) 

[H,D 21 (X)\ = -2(e ip W + e - ik ^)D 21 (X). (5.16) 

The above results justify the interpretation of £> a i(A),C 2a (A) and D 2 i(X) as creation opera- 
tors. B a i(X) and C 2a (X) create single particle excitations, whereas D 2 \{X) creates a bound 



state of two particles. Let us investigate several examples. Since H\G) = 0, (5.14) implies, 
for instance, 

HC 2a (X)\0) = -(2cosMA) + [//2)C 2a (A)|0), (5-17) 

or more generally 

n 

HC 2 a 1 (X 1 ) . . . C 2a „(A n )|0) = - £(2 cos k(Xj) + U/2) C 2ai (Ai) . . . C 2a „(A n )|0). (5.18) 

j'=i 

A similar result holds for states, where the operators i?i a (A) or mixed products of operators 
B a \(X) and C 2a (A) are applied to the vacuum. We have to remember here that B a i(X) and 
C2<x(m) create bounded states, only if p(X) and are real. The restrictions on fc(A) and 
p(X) occurring in ( |5.16| ) have been, k = q — in, p = q + in, where q is real and n > 0. Taking 
these restrictions into account we obtain, for instance, 

HD 21 (X)\0) = -(4cosgcoshK + U)D 21 (X)\0). (5.19) 

The constraint ( |5.7|) implies 



4cos gcosh k = dby 16cos 2 q + U 2 , (5.20) 

where the plus sign has to be taken for U < 0, \q\ < ^, and the minus sign is relevant, if 
U > and ^ < \q\ < ir. This is, of course, in accordance with intuition. 

5.2 Higher Conserved Quantities 



Within the present formalism it is of course more natural to consider the quantum determi- 
nant Det 9 (j4(A)) rather than the Hamiltonian. It can be seen from the commutation relations 
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in Appendix B that arbitrary products of operators B a \(X), C2 a (A) and D2i(X) create eigen- 
states of the quantum determinant Det 9 -A(A). For example, ( [B.26D and Det q A(p,)\0) = |0) 
yield 



Det 9 (A(/i))C 2o (A)|0) 



C 2a (X)\0), 



(5.21) 



pg(X,n)pg(X,jl) v(X)-v(ji) 
where ft is defined by v{p) = v(p) — iU. The eigenvalue of multi-particle states is a product 
of eigenvalues of single particle states. For this reason lnDet q (A(p,)) is a generating function 
of conserved operators, which is possessing an additive spectrum on multi-particle states. 

Commuting conserved operators for the Hubbard model have been constructed by several 
authors p9, 19, 44, 45] either by use of ad hoc methods or by using the expansion (2.8). 



We present the first few known of them in Appendix D. Here we ask for the relation of 
these known conserved operators to the ones generated by ]n~Det q (A(p)). To this end let us 
compare the action on one-particle states. Equation (5.21) implies 



lnDet ? (A(/,))^e-^4|0) = T(p, k) ^ e _ifcj c] CT |0), 

3 j 



where 



T(fi,k) =ln 



k + p(p) k + pijx) 
cos g COS 2 sin£;-sin£;(^) - u7/2 



sm ■ 



k(u) k 
sm — 



sin k — sin k(p) 



(5.22) 



(5.23) 



2 2 

Comparing the expansion of T (//,/;) in terms of v(X)~ 1 with the eigenvalues of the first 
few explicitly known higher conserved quantities H\(= H) , H2, H3, (for details see 
Appendix D| ), we are led to the following conjecture, 



lnBet q {A(p)) 
iU 



+ 



+ 



iU 
iU 



+ 



Lffi + 
i^3- 
iH 4 - 

1 



iU 
2 

- 2UH 3 + 



(iH 2 - UHx 



3U 3 



6U H 



v(p) 



(5.24) 



The coefficients of the asymptotic expansion of In T)et q (A(p)) are linear combinations of 
the formerly known conserved quantities. Particularly, they do not give a complete set of 
conserved operators of the free fermion model in the limit U — > (cf. Appendix D ). It 



21 



remains therefore an open question, if lnDet g (j4(A)) generates a complete set of commuting 
operators for the Hubbard model or not. Another fact which can be concluded from ( [5,24 ) 
is the Yangian invariance of the operators H s (cf . |47]] ) . 

At this point a remark on the construction of Dunkl operators may be in order. Dunkl 
operators are commuting difference (differential) operators, which turned out to be useful in 



the investigation of exactly solvable long-range interacting systems 48, 38]. More precisely, 
a set of difference operators {dj} is a set of Dunkl operators, if there exists a representation 
{Kij} of the symmetric group, such that the operators dj,Kij form a representation of the 
degenerate affine Hecke algebra. Given a representation of the degenerate afflne Hecke alge- 
bra, it is possible to construct a corresponding representation of the Y(su(2)) Yangian out of 



it |48|, [3£| . The strategy developed in ref . [[l^, £38| was successfully applied to the fermionic 
nonlinear Schrodinger model |25|]. Yet it seems to be inappropriate for the Hubbard model 
for the following reason. Let Poj be a permutation operator acting on su(p) spins. Then it 
follows from the defining relations of the degenerate affine Hecke algebra that the transfer 
matrix 



f (u) 



1 + 



ic R 



01 



1 + 



ic B 



On 



d n 



(5.25) 



u — d\ , 

preserves the space of fermionic wave functions. It satisfies the Yang-Baxter relations with 
i?-matrix u + icP and thus generates a representation of Y(su(p)). In case of the fermionic 
nonlinear Schrodinger model the quantum determinant Det q (T(u)) agrees with the quantum 
determinant of a submatrix of the monodromy matrix obtained within the quantum inverse 
scattering approach 123] . It has the eigenvalue 



n 

3=1 



1+ 



IC 



u — k 



(5.26) 



,1 



By way of contrast the eigenvalue of the quantum determinant T>et q (A(^i)) is 

kj + p(jl) 



kj + p{n) n j 
cos — cos ■ 



sinfej — sin k(fi) — iU/2 



sin 



kj — k(fi) . kj — k(jl) sin kj — sin k(/j.) 



sin 



(5.27) 



Hence it seems that the method developed in ref. |^5| has to be modified, if we want to apply 
it to the Hubbard model. 
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5.3 The elements of the monodromy matrix under Yangian transforma- 
tions 

As we have shown above, the submatrix A(X) of the monodromy matrix T(X) generates a 
representation of Y(su(2)). Let us look for the commutators of the remaining elements of 
the monodromy matrix, which can be arranged in submatrices B{\), C(A), D(X), with the 
Yangian generators Q^, n = 0,1] a = x,y,z. Combining ([4.5D and ( |B.8| )-( p3^ ) we end up 
with 

[Q a o,B(\)} = -\a a B{\), (5.28) 

[QJ,B(A)] = smp(X)a a B(X) + ^e abc a b B(X)Q c , (5.29) 

[Qo,C(X)} = l -C{X)~a a , (5.30) 

[Q a x ,C{X)] = -smk(X)C(X)a a + ^e abc C(X)a b Q c , (5.31) 
[Q a ,D(X)] = [Ql,D(\)] = 0. (5.32) 

In the next section we will see that these equations determine the behavior of the eigenstates 
of the Hamiltonian under Yangian transformations. The discussion of the irreducible repre- 
sentations on the subspace of a fixed number of one-particle excitations created by C2 a (A) 
can be done in analogy with ref. [25]. 



6 Construction of n-particle states 
6.1 Scattering States 

We have seen in the preceding section that the repeated action of operators £?i a (A), C2 Q (A) on 
the vacuum produces n-particle eigenstates of the quantum determinant of A(X). For small 
enough n the corresponding wave functions can be worked out by hand. They are of the 
form of Bethe wave functions and are easily understood as scattering states of n-particles. 
For scattering states there is a natural normalization. We have to require the amplitude 
of the incident wave to be unity. In previously studied cases [26, 27, 49 1 it turned out that 



such kind of normalization was obtainable by introducing the operator analog of the reflection 
coefficient of the corresponding classical problem. In other words, the creation operators were 
normalized by multiplying with the inverse of certain generators of conserved quantities. 
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For the Hubbard model we propose the following two pairs of normalized creation oper- 
ators, 

i? a (A)t = i 3 -"e' l Wco S AC 2a (A)Z) 22 (A)- 1 (a = l,2), (6.1) 
£ Q (A)t = p-V^WcosA J B 3 „ Qi i(A)D n (A)- 1 (a = 1,2). (6.2) 

In these formulae a = 1 corresponds to spin-up and a = 2 to spin-down, respectively. The 
numerical prefactors have been obtained by demanding that R a (X) ji and iZ Q (A)^ generate 
normalized one-particle states above the vacuum, 

^(A) f |0) = £ e ~ imfc(A)c U°>> iUA)t|0) = E^ imp(A) 4, Q |0). (6.3) 

m m 

Hereafter we assume that A is chosen in such a way that R a (X)^ and R a (Xy create physical 
states. This means for R a (Xy that k(X) has to be real and for R a (Xy that p(X) has to be 
real. 



It is not difficult to see that eqs. ( 2.29[ )-( p.32 ), which determine the behavior of the 



elements of the monodromy matrix under hermitian conjugation, remain valid in the ther- 
modynamic limit. Hence, they can be used to obtain the conjugated annihilation operators 
corresponding to R a (X) j ' and i? Q (A)^, 

R a (X) = i 2 - a e h W S mX' DniXT'Cij-aiX') (a = 1,2), (6.4) 
R a (X) = i a - 2 e- h ^'hmX' D 22 (X')- 1 B a 2(X') (a = 1,2), (6.5) 

where A' = n/2 — A*. The commutation relations between the normalized operators are easily 
calculated by use of the formulae presented in Appendix B| . Provided A / /i (mod 27r) the 
results are 

R a (X)^R p (^ = -r(A,/i) 75 , Q/3 i? 7 (/x)ti? (5 (A) t , (6.6) 
R a (X)R p (^ = -r^,X) jaj spRy(^Rs(X), (6.7) 
R a (X^R p (p)^ = -r(/i,A) 75 ,^i? 7 ( A i) t i? 5 (A)t, (6.8) 

R a (X)R p (^ = -r(\,ii) y0hS pRy(M?R8W> ( 6 - 9 ) 
i? a (A)t J R /3 ( M )t = -^( M )t J R a (A)t, (6.10) 

R a (X)R p (^ = -Rp(^R a (X). (6.11) 

The operators R a (X), R a (X)^ and R a (X), R a (Xy form a representation of the graded 
Zamolodchikov-Faddeev algebra with 5-matrix r(A, /i). These representations may be identi- 



fied as representations of left and right Zamolodchikov-Faddeev algebra, respectively [50, 26, 
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[51L pL p2| . The grading is such that all operators are odd. In physical terms we may say that 
R a (X)^ and R a (X)^ are creation operators of fermionic quasi-particles. Both quasi-particles 
are charge density waves, because R a (X)^ and R a (X)^ add a particle to the system. Note 
that we cannot have spin density waves over the zero density vacuum. 

It is possible to calculate the S-matrix for two-body scattering processes within the coor- 
dinate Bethe ansatz from two-body phase shifts |5^, |54[] . This approach was applied to the 
Hubbard model at half- filling by Efiler and Korepin |p6f| . The presence of a finite density 
background of particles influences the scattering. Thus the 5-matrix of Efiler and Korepin 
for holon-holon scattering differs from ours by a dressing factor. 

Let us make our above statements about the creation of normalized scattering states more 
precise. We shall present the two-particle states generated by R a (X)^ or R a (X)^ as derived 



from ( 3.11 ) and the commutation relations between the elements of the monodromy matrix. 

fliCAMootio) = E4 T c L T e " infc(A)e " imfc(M) i°)' ( 6 - 12 ) 



R 2 (X)^R 2 (^\0 

E 1 (A)t J R 2 ( M )t|0 

+8(n < m 
R 2 (X)^R 1 (^\0 

+8(n < m 
jRi(A)tAi(/x) t [0 
R 2 (X)^R 2 (^\0 

E 1 (A)t J R 2 ( M )t|0 

+8(n > m 
B 2 (A)t4(/i)t| 

+#(n > m 



E 



C nl C ml e 



J J 



-ink(X) —imk(fj.) 



|0), 



(6.13) 



E C nTS i; 



9{n > m)e -»^) e - im MriJ^!M_ 

v(X) - v(fj,) + iU 



<:■ c '~ J + 9(n < m)e ' c 



-infc(A)g— imk(fi) 

z E C ni c n 
n,m 

-infc(A)g— imk(fi) 



imk(\) — ink(fi) 



-iU 



v(X) — v(fi) + iU 



|0),(6.14) 



Bin > m )e~ ink ^e- imk ^ " (A) 

v(X)-v(fi)+iU 



c C + (H n < iii )e " c 



imk(X) — ink(fj.) 



-iU 



v{X)-v{fi) + iU 



T c mT e -i«pM e -i™P(M)| ) 5 



E 



C il C Ll e 

del 



-inp(X) -imp(/i) 



|0>, 



|0),(6.15) 
(6.16) 

(6.17) 



E nTWni 



6(n < m )e-^(A) e -im P (M) J£LZ?M 

v(X) — v(fi) — iU 



c -\np(X) e -irap(p) + q, > ,-imp{X) ^-inp^) 



iU 



v(X) - v(p) - W 



|0),(6.18) 



E r ' 



.t J 
ni ml 



9{U < m)e -in P (A) e -im pM V(X)- V^) 

v(X) — v(fi) — ill 



-inp(X) -imp(p) 



v ■ ■■<■ -v e "-w/ + 0( n > m ) e "'wv e 



imp(A) „-inp(/i)_ 



iU 



v(X) - v(p) - iU 



|0).(6.19) 
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Note that the two-particle states ( |6,12j )-( 6,15 ) generated by R a (Xy are in-states if k(X) 



< 



k(fi) and out-states if k(X) > k(p,). Moreover, they are normalized in the sense explained 
above. As for the operators R a (X)^ we observe similar things. The two-particle states ( |6.16| )- 
( 6. 19| ) are normalized in-states if p(A) > p(p) and normalized out-states if p(X) < p(n). These 



facts, together with the examples of other integrable models [27|, [4!| lead us to the following 
conjecture: 

Conjecture 1 Provided k(Xj) is real for j = 1, • • • , n, the n-particle state 

i? Ql (Ai) t ---i? Qn (A n ) t |0) (6.20) 

is a normalized in-state if k(X\) < ■ ■ ■ < k(X n ) and a normalized out-state if k(X±) > ■ ■ ■ > 
k(X n ). 

Provided p(fJ>j) is real for j = 1, • • • , n, the n-particle state 

R ai (Vl) j ■■■RaAVnYlO) (6.21) 

is a normalized in-state if p(p\) > ■ ■ ■ > p(^ n ) an d a normalized out-state if p(fJ,i) < ■ ■ ■ < 

The proof of this conjecture seems difficult for general n, since it seems to be unavoidable to 
use the series ( 3.11 ) and the explicit form (|3.13|) of C m (X). 



We have two pairs of normalized one-particle creation operators now, but as in the case of 
B a i(X) and C2 a (A), we do not need to care about both of them in constructing multi-particle 
states. We may use the operator R a (X)> only (or R a (X) only). The reason is the following. 
From fl6.3| ) we deduce that 

fl a (A)t|0) = E a (A)t|0), (6.22) 

where p(X) = k(X). Hence the action of a mixed product of R' and W on the vacuum can 
be expressed in the form ( 6.2C| ) by use of ( 6.1C| ) and ( |6.22| ). In particular, one easily obtains 



-v r a i n(n — l) ~ i ~ i 

J R a „(A n ) t --- J R Q1 (A 1 )t|0) = (-l)-3- L /2 ai (Ai)t.../2 aB (A n )t|0) ) (6.23) 

where p(Xj) = k(Xj). The order of the operators is reversed when written in terms of R&(X) 
instead of w a (X). 



It turns out that the introduction of the prefactors in eqs. (6.1) and (|6.2j) also removes 



the twist fl2.22 ) from the commutators of our redefined creation and annihilation operators 
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with the Yangian generators. Using (|5.2Sj )-(5.32) and some of the formulae in Appendix B, 
we obtain 



<■../•• \ -l%(A)+a§ a; (6.24) 

U 



[Q?,i2a(A)t] = - S mk(\)Rp(\Ya% a + £ abc R^a b a Q^ (6.25) 



[Q a o,Ra(\^] = lM\)^ a p a , (6.26) 
[Q?, i? Q (A)t] = - sinp(A)^(A)t^ Q - ^e abc Rp{\)^ b pa Ql. (6.27) 

These formulae induce an adjoint action of the Yangian on n-particle states |25|. Noting 
that Qo|0) = = Qf|0), we obtain the action of the Yangian on the n = 1 sector as 

Q a R a (X?\Q) = i^ a ^(A)t|0), (6.28) 
Q? J R Q (A)t|0) = -sinA;(A)^ a ^(A)t|0). (6.29) 

Since the action of Qf is — 2 sin k(X) times that of Qq, the representation is called the funda- 



mental representation W\{— 2 sin k(X)) |4(], 41]. 
In the two-particle sector (n = 2) we get 

Qgi? Q (A 1 )t J R (7 (A 2 )t|0) = ^J pa + ^ a a^R p (X^R p (X 2 ^\0), (6.30) 

Q? J R a (Ai) t ii (T (A 2 )t|0) = (-smA;(Ai)^ a V - sin k(X 2 )dp a a a pa + ^£*<7^j 

•i? /3 (Ai) t J R p (A 2 ) t |0). (6.31) 

This representation is a tensor product representation W\[— 2sin&(Ai)) (8> Wi(— 2 sin fe(A 2 )) 
with co-multiplication A defined by 

A(Qq) = Qg®l + l®Qo, (6.32) 
A(Q?) = Q? ® 1 + 1 ® Q? + Ue abc Q b ® Q c . (6.33) 

It is four-dimensional and irreducible, since k(X±) and A;(A 2 ) are real. Due to the Yangian 
invariance of the Hamiltonian, these four states are degenerate. Under the sub-algebra su(2) 
of spins this multiplet is decoupled to su(2)-triplet and su(2)-singlet. We can say that the 
Yangian Y(su(2)) mixes spin multiplets to form a larger multiplet. 

Similarly, the n-particle states i? ai (Ai)' • • • R an (X n y |0) (ay = 1,2) transform under 
Y(su(2)) as tensor product representations W\{— 2 sin k(X\)) <S> ■ • ■ <8> W\(— 2 sin k(X n )). These 



27 



representations are irreducible under the Yangian Y(su(2)), since the quasi-momenta k(Xj) 
are real, but not irreducible under the sub-algebra su(2). In other words, these 2 n states form 
a multiplet under Y(su(2)), while under su(2) they decay into some multiplets according to 
the value of the total spin. Thus su(2) is not sufficient to explain the large degeneracy of the 
system in the thermodynamic limit. 

The irreducibility leads us to the conclusion that we can construct all the ra-particle states 
(|6.20| ) out of the Yangian highest weight state 

R^X^ ■■■R 1 (X n y\0), (6.34) 

by acting with Yangian generators Q^- The wave function of the above state ( |6.34| ) must be 
of plane-wave form, since the on-site interaction never occurs between up-spin particles due 
to the Pauli principle. Therefore, assuming that the state (6.34) is a normalized one (see 



Conjecture 1), we conjecture that the above state ( |6.34 ) is equal to the state 



c\(k(X 1 ))---c\(k(X n ))\0), (6.35) 

where cl(k) = J2j c ] a e ~^ k - Thus we have got a simple method for constructing multi-particle 
scattering states. They are obtained out of the plane- wave state ( |6.35 ) by using the Yangian 



generators ( f4.12|) and ( [4.13] ) . We have already encountered such kind of situation in the case 
of the repulsive <5-function fermi gas p! 



One can similarly discuss Yangian representations of multi-particle states constructed by 
use of W a (X). The alert reader will have noticed the different signs in front of U in eqs. ( |6.25D 
and ( 6.27[ ), which lead to different definitions of the co-multiplication (cf. ( |6.32 ), ( 6.33| )): 



A'(Qg) = Qg®l + l®Qg, (6.36) 
A'(Q?) = Qi ® 1 + 1 ® Qi - Ue abc Ql ® Q c . (6.37) 



But this does not cause any contradiction. The order of the quasi-momenta in (|6~23| ) is 



reversed in the multi-particle states expressed by R a (X)^ compared to those expressed by 
R a (X)^ . This corresponds to the reversed order of the tensor product (g> in the definition of 
the co- multiplication, which compensates the different sign in front of U in ( |6.33| ) and (3.37). 



The various parameters U, A, h, v, k and p which we used so far are connected through 



the formulae (|2.5|) and (4.3). Thus only two of them are independent. As a test of consistency 



of the results in this section let us consider the free fermion limit U — > 0. This limit is most 
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conveniently taken for fixed v, for if we fix v = v(X) and v = v(fi) in eqs. ( |6.12| )-( |6.19| ) and 
let U approach 0, we see that the products of operators R a (\y and R a (\y act like products 
of creation operators of Bloch states c4(/co) on the vacuum. Here ko = po is determined by 
the corresponding limit in eq. ( |4.3| ), 

v 

sin ko = — — . (6.38) 



A and h are now dependent variables. Considering (2.5) and (4.3) for fixed v and small U we 
find the following solutions 

11 ( v 2 \~^ 

icotA = 1 + T\ 1 ~t) +°( [/2 )' ( 6 - 39 ) 
e 2h = i(l- V -pj 2 - V - + 0(U 2 ). (6.40) 

Using these equations and some standard trigonometric identities we can express all the 
functions of h and A, which enter the definition of C m (X) (cf. fl3.13|) ) in terms of v and U. 
Note that ( 6.3S| ) and ( |6.4[1| ) are not the only possible solution of ( |2.5|) and ( |4.3| ) for fixed 



v 



and small U. We choose the branches such that \\mjj^Q C m {X)\ v = For small U the odd 
elements of C m (X) — I4 are of the order of and the even elements are of the order of U. 



Thus only the first sum on the rhs of (3.11) contributes in order Uz to the odd elements of 
T(A) — I4, and we obtain 

C 2a (\) = r-*^J24n«e- imko + 0(U 3 2), (6.41) 



in 



B 3 - a ,i(X) = i^-^-^c^e-^ +0(UI), (6.42) 
cosA m 

where e ±h / cos A = 0(U^). Since Dpp(\) = 1+0 {U) (f3 = 1, 2), it follows from the definitions 
flOP and Q that 

lim R a {X) ] \v = hm R a {X) ] \v = cl(k ). (6.43) 

The corresponding formulae for R a (\) and R a (X) are true by hermitian conjugation. Eqs. 
- Q6.11 ) turn into the usual anticommutators between fermi operators, since 



limr(A,/i)| v - = P. (6.44) 

We see that we may interpret the Zamolodchikov-Faddeev algebra as a deformation of the 
anticommutators between fermi operators with deformation parameter U. 
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6.2 Bound States 

One of the delicate points in Bethe ansatz calculations is the question of completeness. The 
completeness of the coordinate Bethe ansatz for the Hubbard model under periodic boundary 
conditions was discussed by Efiler, Korepin and Schoutens |33|. They showed that the Bethe 
wave functions are highest weight with respect to the su(2)©su(2) symmetry generated by S a 



and rf [33]. Hence, each solution of the Bethe ansatz equations correspond to a su(2)©su(2) 
multiplet. They counted the number of Bethe ansatz solutions assuming Takahashi's string 
hypothesis (3^] to be valid and multiplied by the multiplicities of the su(2)©su(2) multiplets. 
The resulting number is equal to the dimension of the Hilbert space. Note however that, 
although a matter of common belief now, the string hypothesis is waiting for a proof since 
25 years. 

How to pose the question of completeness in our infinite chain formalism? As we have seen 
in the preceding section the operators R a (\y and R a (X)^ create single electrons in scattering 
states. The operator L>2i(A) creates a bound pair of electrons. As we have learned in section 
5, there are no operators that create more than two particles among the elements of the 
monodromy matrix. As we will see, however, the string hypothesis suggests the existence of 
bound states of pairs. How to define the corresponding bound state operators? 

To obtain a guess, let us recall the string hypothesis. In the following we will denote the 
spin rapidities of the coordinate Bethe ansatz by Aj and the momenta by kj (for details cf. 
[HH). According to the string hypothesis there are two types of string solutions of the Bethe 
ansatz equations. 

1. (A-string) m Aj's form a string configuration, in which the real parts of the A^-'s are 
the same while the imaginary parts are arranged at equal spacing of ill /2. The center 
of the string should be real. 

2. (k- A-string) 2m /-q's and m A^'s form a string configuration. The values of ki's and 



Aj's are 



k\ = 7T — arcsin(A' + imC//4), 
k 2 = arcsin(A' + \(m - 2)11/4), 
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&2m-2 = arcsin(A' — i(m — 2)t//4), 

k2m-l = 7T - &2ro-2 

^2m = vr — arcsin(A — imU/A), 

Aj = A' + i(m + 1 - 2j)U/4, A' real, j = 1, 2, • • • , m. 

These solutions should be exact in the thermodynamic limit. Since we are dealing with the 
zero density vacuum, there should be no spin excitations, and we do not have to consider the 
A-string here. 

We wish to obtain an operator which creates a 2m-A;-A-string (for short we shall simply 
call it as "2m-string" ) . To begin with, we shall deal with the 2-string (m = 1 case). There are 
2 particles involved, one with spin up and the other one with spin down. The wave function 
of the 2-string state is of the form 

E c lAl { 6{m ^ n)e-^ mk+n ^ + 9{m < n )e^ nk+m ^} |0), (6.45) 

m,n 

where sin k — sin p = iU/2 and k+p real. Therefore, by comparison with eq. ( |5.6| ) , the 2-string 
state is proportional to L>2i(A)|0) with an appropriate choice of A. It can be easily seen by 
explicit use of ( |3~Tl| ) and fTiiP that it is also proportional to C 22 (A 7 )^! (A") |0) , if A' and A" 



satisfy the following conditions; 

p(X') = vr - k{\") mod 2tt, (6.46) 

k(\") = p(X) mod 2vr, (6.47) 

jfe(A') = fc(A) mod 2vr. (6.48) 

These are three conditions for three parameters A, A', A", which at first sight seems to violate 
the arbitrariness of A. Yet there is a redundancy in these equations. ( |6.46 ) and ( 6.47] ) imply 



p(A') = vr - p(X) mod 2vr, (6.49) 

which is compatible with (|6.48| ) by taking into account the constraint ( |5.7j ). Thus we have 
obtained two possible 2-string creation operators, which are connected with each other by 

ie h(X')+h(\") CQS A // CQS 2 A / sin y 



D - (A>|0) " i^A 

1 _ e i(p(A')-*(A"))Vl - e im")-k(X))) 

J -C 2 2(X')C 21 (X")\0). (6.50) 



1 _ e i(p(A')-fc(A')) 
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A, A' and A" in this equation have to satisfy (|6.46|) - (|6.48| ). Note that it follows from ( p.5| ) 
that 

C 22 (A')C 21 (A") = -C 2 i(A')C 2 2(A"). (6.51) 

Let us proceed with the general 2m-string states. We conjecture that the creation operator 
of a 2m-fc-A-string can be expressed as 

C^ m \\i, • • • , A 2m ) = C 22 (Ai)C 2 i(A 2 )C 22 (A3)C 2 i(A4) • • • C 22 (A 2m _i)C 2 i(A 2m ), (6.52) 



where 



k(\2s) +p(A 2s _i) = vr (mod 2ir), 

. , ,. . . , . iU(s — 1) , 
sinfc(A 2s -i) = smfc(Ai) H , (s 



1, • • • ,m). 



(6.53) 
(6.54) 



Following previous works [56, [)7], 58, 59 1 we shall call this operator bound-state operator. The 
expression on the rhs is a formal one and should be interpreted as a "composite operator" (see 



Appendix C). One can easily verify that the functions sin/c(Aj) form the same configuration 
as in the /c-A-string, if their center 



C= — > smk(\i) = sinfe(Ai) + 



2m 



i=l 



4 



(6.55) 



is real. 



We can normalize the bound-state operator by a method similar to that in the case of 
the scattering states created by R a (\)< . Let 



£>g m) (Ai, • • • , A 2m ) = J D 22 (A 1 ) J D 22 (A 2 ) J D 22 (A 3 )Z? 22 (A 4 ) • • • I? 22 (A 2m _i) J D 22 (A 2m ) 
We define a normalized bound-state operator as 

R^ 2m \\i, • • • , A 2m )^ = C2 2m \\i, • • • , A 2m )-D 2 2 m (Ai, • • • , A 2m ) 1 . 



(6.56) 



(6.57) 



Similar definitions of bound state operators have appeared before in the literature. They 
have been applied to the XXZ-chain [58, ^] and to the attractive ^-function gas [56, [60| [. 
Note that it is not a priori clear, if we can apply the commutation rules for the elements of 
the monodromy matrix to obtain the commutators of 2m-string operators (cf. Appendix Q ). 
However, if we assume we can, we get the following reasonable result, 

R i2m \Xi^R {2n Hji^ = (~ r l + ( n + m ) iU / A C~r] + \n-m\iU/A 
3 £ — r\ — (n + m)iU '/4 £ — r\ — \n — m\iU/4 
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min(m,n) — 1 _ „ , _ s. Tr/ .-,2 

C - V + {n + m - 2s)i£//4 n 



n 



s=l 

^ 2m )(A,)ti? a (^)t 



C-r/-(n + m- 2s)H7/4 



J R^ 2T1 >(^ J )t J R( 2m )(A i ) t , (6.58) 



£ — sin k{ji) + \Um/A 



R a (p)lR( 2m \\. 



t 

1 1 i 



(6.59) 



C — sin k(fi) — i[/m/4 

where C is the center of the 2m-string and r/ is the center of the 2n-string. The factor on 
the rhs of ( |6.58j ) is the ^-matrix between a 2m-string and a 2n-string, and that on the rhs of 
( 6.59| ) is the S'-matrix between a 2m-string and a particle. ( |6.58| ) is of the same form as the 



5-matrix for the scattering of bound states of magnons in the XXX-chain [ 59 1 . 
As for the transformation under the Yangian Y(su(2)) we can easily show that 

[Q%,R( 2m \\rf} = = [QlR^HXi?], (6.60) 



which follows from eq. dC.6|) and from the commutativity of and -D22(Aj). From ( |6.60| ) 
the action of the Yangian Y(su(2)) on a 2m-string state can be derived as 



Q a Q R^ (A»)t |0) = 0, QfR^ (Ai) f |0) = 0, 



(6.61) 



i.e. the 2m-string state is singlet under Y(su(2)). 



7 Concluding Remarks and Discussion 

We have developed the QISM for the Hubbard model on the infinite interval with respect 



to the zero density vacuum. The i?-matrix ( 3.27 ) thus obtained is greatly simplified in 
comparison with the i?-matrix of the finite periodic model. In particular, it reveals a hidden 
rational structure, which arises from a certain combination of the functions pi in eq. 



This structure was discovered earlier by Ramos and Martins [22| as part of the exchange 
relation for the Hubbard model on the finite interval. Along with the simplified i?-matrix we 
obtained the asymptotic expansion ( |4.5| ) of the submatrix A(X) of the monodromy matrix, 
which naturally provides a representation of Y(su(2)) and generates an infinite series of 
mutually commuting Yangian invariant operators, which is including the Hamiltonian. We 
thus clarified the origin of the Yangian symmetry of the Hubbard model. 

We constructed creation and annihilation operators of elementary excitations. There are 
two types of excitations over the zero density vacuum, one-particle excitations and 2m-string 
excitations. We showed that they are distinguished by their behavior under the action of 
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the Yangian. All 2m-string excitations are Yangian singlets, whereas the n-fold one-particle 
excitations transform like n-fold tensor products of the fundamental Yangian representation. 

The interaction of two elementary excitations is described by their S-matrix, which is 
given by the commutation relations of the corresponding operators. We calculated these 
commutation relations and thus the S'-matrix. The 2m-string excitations are Yangian sin- 
glet, which means that they have no internal structure. Therefore their mutual interaction 
produces only phase shifts, their S'-matrix (6.58) is a scalar factor. In contrast, the scatter- 
ing of one-particle excitations involves the spin. The S-matrix is a 4 x 4-matrix. Up to a 
factor of minus one, which accounts for the fermionic nature of the excitations, it agrees with 
the rational submatrix r(A,/i) of the i?-matrix in the thermodynamic limit. We obtained 
two alternative pairs of one-particle creation operators. The corresponding S-matrices are 
reciprocal to one-another, which means that, if one of the pairs creates in-states, the other 
one creates out-states. Our creation and annihilation operators combine into a left and right 
representation of the Zamolodchikov-Faddeev algebra. Although it should be obvious, we 
would like to emphasize that our one-particle operators in equation (6.6)-(6.11) are not ab- 
stract quantities, but explicitly given in terms of the elementary creation and annihilation 
operators of Wannier states. In the limit of vanishing coupling (U — > 0) our one-particle 
operators turn into creation and annihilation operators of Bloch states. The Zamolodchikov- 
Faddeev algebra can thus be understood as a deformation of the anticommutators between 
fermi operators with deformation parameter U. 



There is a known standard method for the solution of the "quantum inverse problem" 27], 
which consists in expressing the fermi operators for electrons in Wannier states Cj a , cj CT in 
terms of quasi- particle operators R a (X), R a (X)^ and R^ 2m \\), i?( 2m )(A)t. In analogy to the 
classical inverse scattering method, we have to derive a quantum Gelfand-Levitan equation. 
To this end we should investigate the analytic properties of the monodromy matrix. Due to 
the complicated structure of the i?-matrix, this may be a difficult task. We leave it for future 
work. The solution of the inverse problem will enable us to calculate Green's function with 



respect to the considered vacuum, as it was done before for other integrable models |3C], 58, 5S]. 

So far our work has been limited to the four cases of uncorrelated vacua, for which the up- 
spin or down-spin orbitals are either completely filled or vacant. Unfortunately, these cases 
are not of particular interest from the point of view of condensed matter physics. Correlation 
functions are rather trivial and may be obtained directly within the coordinate Bethe ansatz 
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solution of Lieb and Wu. Hence there is a clear need to develop a method for the algebraic 
construction of excitations over a correlated vacuum. In a first step we had to renormalize 
T mn {X) with respect to this vacuum. The properties of the vacuum enter the formalism, they 
are not results of it, and therefore have to be obtained by independent means. Assume we 
had done this step. Then the next question would be, if the vacuum is a Fock vacuum with 
respect to some of the entries of the renormalized monodromy matrix. This is, of course, 
not clear a priori, but is a necessary requirement for an algebraic construction of eigenstates 
based on commutation relations. 

The development of a version of the QISM for interacting vacua is certainly a hard and 
challenging task, and it may be more appropriate to start with a more simple model than 
the Hubbard model. It seems however, that such kind of method would finally lead to a 
complete understanding of the Hubbard model, too. We would like to consider our work as 
a first necessary step towards this goal. 

For a further investigation of the Hubbard model by algebraic means the physically inter- 
esting half filled case seems to be most appropriate. At half filling the 2m-strings disappear 
from the spectrum 46]. There are only two pairs of independent one-particle excitations, 



which have been described as "spinons" and "holons" HqJ. In analogy to the results in sec- 
tion 6.1 we expect the space of states for a given set of rapidities to be spanned by the 
action of Y(su(2))©Y(su(2)) Yangian generators on Yangian highest weight states. We fur- 
ther know that the half filled ground state (at zero magnetic field) is singlet with respect 
to Y(su(2))©Y(su(2)) f42|| . This suggests that form factors may be calculated by purely 
algebraic means Q. 
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Appendix A Singular Terms in the Infinite-chain Formalism 

In this appendix we calculate the limits 

{U + {\,^Y l ) a(3a5 = lim (E/ n (A,/i) -1 ) a/ 3 7 ,5, (A.l) 
U-(\,fi) a/ 3js = lim U n (X,fjL) a 8 7 s, (A.2) 



which determine the elements of T(A) (g> T{p) via equation ( 3.24| ). U n (X, fi) 1 and U n (X, fi) 



are defined by eq.( |3.22 ), where V(\) is given according to fl3.12| ). V^(X, p) is easily obtained 



by direct calculation. Its diagonal consists of the elements of V(X) ® s V(p,). Due to normal 
ordering there appear additional non-vanishing off-diagonal elements, 

^ (2) (A,^) 12i21 = F (2) (A,/i)i 3 ,3i = -i sin A sin /i, (A.3) 
^ (2) (A,^)i4,23 = -^ (2) (A,^)i4,32 = -isinA cos M , (A.4) 

F (2) (A, /U)24,42 = ^ (2) (A, ^34,43 = i COS A COS fl, (A.5) 

^ (2) (A, m)23,41 = (A, /i) 3 2,4i = -i cos A sin M , (A.6) 

U (2) (A, ^14,41 = -e h M +h M . (A.7) 

Note that V^(X,fi) is upper triangular. Since the diagonals of y( 2 )(A,^) and V(X) (8> s V"(/i) 
agree, V^ 2 \X,fi) can be diagonalized by an upper triangular matrix U(X,p) whose diagonal 
elements are all unity, and 

V( 2 \X^) = U(X,ri(V(X)® s V( f ,))U(X^)- 1 . (A.8) 

It turns out that the non-vanishing off-diagonal elements of U(X, n) are simple rational func- 
tions of the Boltzmann weights pj = pj(X,p). They are obtained as 

U(X, M)l2,21 = U(X, /i)i 3 ,31 = — , (A.9) 

Pw 

U(X, M)l4,23 = -U(X, M)l4,32 = — — , (A.10) 

P3 - Pi 

U(X, /i) 24 ,42 = U(X, p) UA3 = ^, (A.ll) 

P9 

U(X, ^23,41 = ~U(X, /i) 32 ,41 = (A.12) 

P5 - Pi 

U{X,ii) UA1 = ^^. (A.13) 

P5 - Pi 

Before proceeding further let us introduce some shorthand notation. Instead of /(A), g{n) we 
will write /, g. The bar means that the argument of the function if p. Using this convention 
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and equation (A. 8) we obtain 

U n (\,v) = U(\,/j,)(V~ n <g> s V~ n )U(X,fi)~ 1 (V n (g> s V n ), (A.14) 

Unix^y 1 = (y~ n <g> 8 v~ n )u(x,fjL)(v n ® s v n )u(x, p)" 1 , (A.15) 

Comparing ( A.14 ), ( |A,15 ) and ( A.l ), ( |A.2| ) we find that we have to calculate the following 
limits 



It 



if 



If 



ii m ]El{\ _jn(p-py 

7W±0O p 1Q 



1/02 

-0 

lim — — 

n— >±oo pg 



lim 
lim 



±00 - p 3 
ipe 



n->±oo p 4 — p 5 
P3 



(1 -e in(fe ~ fc) ), 

/]_ _ e in(p-fe)\ 
(1 _ e i«( fc -p)^ 



(A.16) 
(A.17) 
(A.18) 
(A.19) 



lim 



+ 



3 in(p-fc) 



00 1 >0l — /03 (Pi - P3)(P4 - P5) 



+ 



P5 j n (fc-|-p— fc_ p ) 



P4 - P5 



(A.20) 



lim 

n— ►— 00 



P5 



+ 



2p 



jn(k-p) 



P3 

H e 



(Pi - P3)(P4 ~ Ph) 
in(fc+p— fc— p) 



(A.21) 

Pi - P3 

These limits exist in the sense of generalized functions. It turns out that they all can be 
reduced to the following formula 

1 _ e i(p-p> 1 



lim 



(A.22) 



n^ioo e [ P — e J P e ip — e^ p±i£ ^ ' 
where e is an infinitesimal positive number. The limits exist in the domain ±Q(p — p) > 
A proof of the above formula can be obtained in two steps. First show by acting on a test 
function that 



lim p.v.—r 

i. — >-\-r>o ci l \ 



n-*±oo * e Kp-p) — 1 

where 5 is the periodic <5-function 



±w8(p-p), 



6(p) = 5(p — 2nir). 

n=—oo 

Then use a periodic version of the Plemelj formula (cf. page 501 of ref. 

., , . , = p. v.- — ± tt6(p), 

1 _ e i(p±ie) 1 1 - e'P 



(A.23) 



(A.24) 



(A.25) 



37 



to obtain flA~2l ). 

Using ( JA.22 ) we find, for example, 



[ e -h-h j — e Kp-p) n 

It 1 = lim — - : — 

sin A sin \i n— >±oo e ip — e ip 

ie -h-h 1 



sin A sin [i e^ - e<J>± ie ) ' 



(A.26) 



and similarly 



■ h+h i 

l± = 1 (A.27) 

2 sin A sin ji e ife - e ^ k±i£ ^ ' 

:„-h+h i 

l f = -^T-- •£ (A. 28) 

sin A sin \x e lk - e l( ~P ±1£ > 

[ e h-h 2 

if = — — ^ ttt-t- (A.29) 

4 sin A sin ^ e 1 ? - e l ( fc±1£ ) 

Note that the prefactors on the rhs of the above formulae when expressed in terms of k, k, p, 
p are all regular at the singularities of the second factors. The remaining limits are obtained 
as 

/+ (l + e^Xl+e^+P)) [ j* + j k 2 1 

5 e ij5 _ gife | e i(k+P) _ e i(fc+p+ie) gifc _ e i(p+ie) J ' ^ ' > 

(1 + + eg+g) / e 1 ^ + e^ 2 1 

6 e ifc _ e ip I e i(S+p) _ e i(fc+p-ie) e ip — gi(fe-ie) J ' \ • ) 

The last two equations require consideration of the consistency of the occurring singularities. 

In particular, in case of ( A.30| ) the following two conditions must be satisfied, (i) ^s(jp) > 
Q(k) if k + p and k + p are real, (ii) Q(k) > Q(p), if k and p are real. These two conditions 
have to be consistent with the constraint (|5.7|). 



Consider the first condition. We can use equation (p.q) to obtain 

9(p — fc) = — arcsinh f- — r-r^ — arcsinh ( —= r-r^r^ • (A. 32) 

Hence ^s{p — k) is positive for positive J7, if ^ < < w, and positive for negative U, if 

< § ■ Note that the same restrictions on the parameters were obtained in section 5.1 as 
conditions on £>2i(A) to create a bound state. 

Consider the second condition above. The constraint ( |5.7| ) has two branches of solutions 
for p as a function of fc, which are fixed by choosing sgn(9(j>)). We may therefore choose 
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> > and (ii) will be satisfied. Moreover, this choice of branch ensures the 

prefactor (e 1 ^ — e lfc ) _1 in QA,3C| ) to be nonsingular as k approaches p on the real axis. 

Equation flOl] ) may be discussed in a similar way. Now the two conditions become (i) 
Q(k) < Q(p) if k + p and k + p are real, and (ii) ^(p) > ^(k), if p and k are real. Condition 

(i) implies the same restriction as above, ^ < \^-^\ < vr for U > 0, and |^^| < \ for f7 < 0. 

(ii) can again be satisfied by an appropriate choice of branch of p as a function of k. 

We are now in a position to consider the weak limits 

(T(A) ® s f{n)) a p aS = (L 7 +(A,/i)~ 1 ) Q/3ie(/ ,T (2) (A,/i) e ^ w C/_(A,/i), ?Pi7( 5. (A.33) 

We have to check, whether the functions on the rhs of this equation exist on a common domain. 
Clearly, this has to be done equation by equation. We find the following 16 combinations of 
functions Z- : 

OlJ^l,}; Ol" ) ^3 }' {hil£ih}' 03^5")^}i Ol '^4 ^6 }> 
O3 ' ^5 > ^4 ' ^6 }> {^3 > ^5 '^2 }' 04>^l}> O4 > ^3 }> 
{^2'^ll' O2 ' ^3 }> {^'^4'^}' 0<jt> ^2 }' {^2 '^4'^}' {^2 » ^2 }' 

It turns out that each of these combinations is compatible. However, they are not compatible 
all together. 



Now we should solve ( A.33| ) for (A, /i) Q /3 i7 5 and insert the result into ( |3.20 ) to obtain 



the commutation relations between the elements of the monodromy matrix. We should do 
that equation by equation and should take care of the compatibility of the domains of the 
occurring functions. This would be a cumbersome task and would, moreover, obstruct the 
algebraic structure of our problem. Therefore we leave mathematical rigor at this point and 
proceed more formally. 



First note that there occur products of the form ifli on the rhs of ( A.33| ). These products 



are not well defined, since the regularization requires two limits, and, when acting on a test 
function, the result will depend on the order of these limits. "This indicates the highly 
singular operator character''^ of some of the elements of the monodromy matrix. In the 



following we will therefore exclude the singular points of the functions if (A. 26) - (A. 31). 



We will assume that k ^ k for k, k real etc. Then we can omit the regularizations ±ie in if, 
*We are citing Sklyanin j28[. 
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and 

U+(\,ii)- 1 = U(\,n)- 1 , U-(\, f j) = U(\, f j). (A.34) 
We insert this result into (|A.33| ), treat ([A.33| ) as a matrix equation and solve for T( 2 )(A, /i), 

f®(\, M ) = U(X, /i)(T(A) ® s T(jj,))U(X, n)- x . (A.35) 

Now ( |3.20| ) implies ( p. 25 ), and the i2-matrix ( p. 27 ) is obtained by using ( |A.34 ) in eq.( 3.26| ). 



Appendix B List of Commutation Rules 

Appendix B.l Elementary Commutators 

In this appendix we provide a complete list of the commutation rules encoded in the exchange 
relation ( p. 25 ) in terms of the submatrices A(X), ■ ■ ■ ,D(X) of the monodromy matrix T(A). 



As mentioned in section ||| these submatrices generate sub-algebras of ( 3,25 ), 



r(A lM )(A(A)®A(/i)) = (A0i)®A(A))r(A,M), (B.l) 
a(\,ti)(D{\)®D(ji)) = {D(»)®D(X))s(X,fi), (B.2) 



^^ r{X ^\B 2a (X))®[B 2a M) 

(B la (fi)\ (B la (X)\ 
\B 2a (^)^[B 2a (X))' 

P4(A,m) (x (B 11 (X)\ (B 12 (fx)\ 

Pl (A, / ,)-p 3 (A^) n ,N \B 2l (X)J ® \B 22 {n)) 

= ( BM ) ft ^ ll(A) ^ 
~U 22 (/x)PU 21 (A)j' 



(B.3) 



(B.4) 



(C al (A),C o2 (A))(2)(C ol (//),C a2 (//)) 

= (c al (^,c a2 (^)^(c al (A),a 

(Cn(A), Ci 2 (A)) ft (C 2 i(/i), C 22 (/i)) 



(C al (p),C a2 ^))^(C al (X),C\ 2 (X))^^r(X,^), (B.5) 



(g 2 i(/x),C 22 (M))®(C , ii(A) i Ci2(A)) n P4 $ X,fi) n M X,?). (B.6) 

Pi(A,/i) - p 3 (A,M) 



The matrix s(A,/x) in (|B.2|) is defined as 



S (A, M ) = diag ( 1, 1 ) V. (B.7) 

Pa- Pb Pi 
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The commutators of B(X),C(X) and D(X) with the submatrix A(/j,) are given as 

l/,Jo(A ^ (C n (A),C 12 (A))®^(/x) = (^(/x)®(C 11 (A),C 12 (A))r(A,/x), 



P4(A,//) 
Pq(A,/x) 



(C 21 (A),C 22 (A))®A(/x) = (A(/i)®(C2i(A),C 22 (A))r(A,M). 



[£»n(A),^( M )] = [£> 22 (A),AM] = 0, 

£>2i A A (//) = — - — - — -A{n)D 2 i{\). 

Pi(A,m) P4(A,m) 

The commutators between the entries of B(X) and C(/x) are 

B al {X)CM = Pio(A,M) 2 Cl6(/ , )Bol(A)> 

Pl(A//)P4(A,//J 

£ a i(A)C 26 (/z) = ^^(^.(A), 
P9(A,m) 

5 a2 (A)C 16 (^) = ^Mc 16 ( Ai )B a2 (A), 
5 a2 (A)C 26 (/.) = - ^fc^ C 26 (p)^ a2 (A). 
Finally, there are the following commutators between B(\), C(A) and the submatrix D(A) 

B al (\)D(») i 1 5 ) £>(p) ( 1 ° ) 

Pl ~ p3 



PlO 





pi 







P9 


£10 





Pi 





P9 


(: 







pi 
Pi— PS 



B a2 (\)D(») = i y « P1 J y P4 P5 1 J S «2(A), 

C la (A)£>(/z) = -if 1 ' I £>(//) f^ ° |C la (A), 



P4-P5 Q 

C 2a (A)D(/z) = -il ^ !!)(/*) I <*> " |C 2a (A). 



pi 


I) 


Pi-Ph 







Pi 





PlO 





_Pi 


Pi 


P4 





PlO 





_P9 


Pi 
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Appendix B.2 Commutators with the quantum determinant 



Eqs. ( |B.8D -( [B?L4 ) imply the following commutators of the quantum determinant Det ? (A(/i)) 



with the remaining entries of the monodromy matrix 

rw ( a( \\n tw pj{X,p)p 4 {X,p) v(X)-v(p) m^w,, u /R ^\ 

Det q {A(p))B al {X) = — r B al (X)Det q {A (p )), B.23 

Plo(A,M)Pio(A,/Lt) f(A) - v{fi) + 1(7 

Det 9 (A(/.))£ a2 (A) = - pi(A|Ai)pi(A>A)T , (A) _ t;(M)+i y ga2(A)Det g (^(/i)), (B.24) 
r> f ^ \\^ ^(A, p)p w (X, p) v{X) - v{p) + ill 

Uet q [A{fj,))Ci a {X) = — t-t — C la {X)T)et q {A(n)), (B.25) 

/3 4 (A,/i)p 4 (A,/i) v(A)-u(/i) 

De,,(A(M))C 2o (A) = _ ^(A.M) W (A,W «M - -M +ig Ca . WDeMit0l)) . (B.26) 

p 9 (X,p)p 9 {X,p) v{X)-v{p) 



D 12 (X)Bet q (A(p)), (B.27) 



Det,(A(/i))Di2(A) = 

p 9 (A, p)p 9 {X, /2) pip(A, ^)pio(A, /2) 

pi(A,/i)pi(A,/i) p A {X,p)p A {X,p) 
Bet q (A{p))D 21 {X) = 

p 1 (X,p)pi(X,jl) pa{X,p)pa{X,p) t\\n.(A(\\ 

77 s 7T-^T 77 \ 7\ 7V-^2l(A)Detg(j4(/J,)), B.28) 

p 9 (A, /x)p 9 (A, //) pi (A, //)pio(A, m) 
[Det (? (^(/x)), J D 11 (A)] = = [Det g (^(/x)), J D 22 (A)] ) (B.29) 

where /2 is defined by v (p) = v(p) — ill. 

Appendix C Composite operators 

For the construction of bound state operators corresponding to the 2m-string in section 
6.2, we have to introduce composite operators, which are formal products of entries of the 
monodromy matrix. Our definition of the 2-string creation operator, for instance, is 

i?( 2 )(Ai, A 2 )t = C 22 (Ai)C 2 i(A 2 )Z) 22 (A 2 )- 1 Z) 22 (Ai)- 1 , (C.l) 

where 

p(Ai) + fc(A 2 ) = 7T (mod 2tt). (C.2) 



Because of the constraint (5.7), this implies 



frinfc(Ai) - sin£;(A 2 ) = l —. (C.3) 
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Thus k\ and fc 2 cannot both be real, and it follows from (5.4) and (5.5) that C22(Ai) and 



C21XA2) in the products on the rhs of (C.l) cannot both create bounded states. They 



are not simultaneously physical operators. Despite this fact, we can provide the product 
C22(Ai)C2i(A2) with physical meaning, when sinfc(Ai) — iU/A is real. Eq. ( A.33| ) implies 



C 22 (Ai)C 2 i(A2) = l^ 2 \X 1 ,X 2 ) 44,14 + 44,32- (C.4) 

The matrix elements 7~( 2 \\i, A2) a/ 3, 7 5 are well-defined through a series representation similar 
to ( 3.1 1| ) . We may take ( |C.4D as a definition of the composite operator on the left hand side. 



Its domain is the domain of Z 3 . It is determined by the condition Q(p(Ai)) < Q(k(X2)), which 
leads to the same discussion as above eq. (0) or below eq. (A. 31). 



In principle, we could iterate the renormalization procedure explained in section 3 and 
in Appendix A . We could define £m\Xi, • • • , A m ) as the /c-fold graded tensor product of L- 



matrices at site m and could introduce its expectation value V( k \X\, • • • , X m ) which governs 
the renormalization of the fe-fold tensor product of monodromy matrices Tmn(Xi, • • • , A^). We 
would obtain a renormalized tensor product T^(X\, • • • , A&) and the commutation relations 
between the entries of T( fc )(Ai, • ■ ■ , X k ) and TW(// X , • • • , (let f^\X) = f{X)). We guess 
that such kind of procedure would solve the completeness problem in a satisfactory way. 
Unfortunately, it seems to be practically impossible to do these calculation, because of the 
increasing dimension of the involved matrices. 

A composite operator is not a mere product of the original operators. So it is not obvious 
whether or not the commutation rules between composite operators follow from iterating 
the commutation rules of its factors, as obtained from ( |3.25 ). Nevertheless, we assumed so 



and investigated some of the consequences of this assumption. This way we obtained the 
/S-matrices ( 6.5§| ) and ( |6.59| ) for composite operators, which look very reasonable. Another 



consequence of such kind of formal procedure is that all 2m-string states are Yangian singlet. 
In case of the two-string this can be seen as follows. Using (|B.5|) with P4,{X2-,X\) 7^ 0, 
pi(X 2 , Ai) = and r(A 2 , Ai) = (1 + V)/2, we get 

C 2Q (A 1 )C 2 ^(A 2 ) = -C 2/3 (Ai)C 2Q (A 2 ). (C.5) 

From this equation one can derive the following formula 

[Q8, C 22 (Ai)C 2 i(A 2 )] = = [QJ, C 22 (A 1 )C 21 (A 2 )], (C.6) 

which indicates that the composite operator C22(Ai)C2i(A2) creates a Yangian invariant pair 
of up- and down-spin particles. 
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Appendix D Conserved Quantities of the Hubbard Model 



Explicit expressions for higher conserved quantities have been derived by several authors [44, 
45]. In contrast to the case of the nonlinear Schrodinger model these quantities are not in one 
to one correspondence to the conserved quantities of the free fermion model. Let us define 
creation operators of Bloch states as 



(D.l) 



In terms of these operators we can define two sequences of operators 
F n = ~j dk sin (n (k + cl(k)c a (k), 
F n = -- J dk cos (n (k + 4(k)c a {k). 



(D.2) 
(D.3) 



All of these operators are mutually commuting, and we have H(U = 0) = iq. Using the 



results of ref. [45] we obtain the following higher conserved quantities. 



H l = H, 

3,v 



+ W X)[(cJ+i,<r^> ~ c 3,* c 3 J rh*)( n 3-° + n i+ l -° ~ 

r 3 

+ ( C }+2,cr C j + l,cr — c] + lj(T C :) ' + 2,(T)(c] + 1) _ cr Cj ) _ cr — _ a Cj + \- a ) 
+ 2 ( C ] + l,cr C j,cr — C ] ia c j+l,cr)( C ]+l,-a C j,-tT — C^Cj+^-o-) 



(D.4) 



(D.5) 



5) ( nj '~ CT 



7.-°" 



i\ 3 

8 

1\ 1 
+ - 



7 

iC/ Z][( C i+3,<r C J> - C jf,ir C J+3,<T) ( n j,-cr + n j+ l- a + n j+2 - a + ™j+3 -a ~ 2) 
+ (4+2,<7 C i,<x + c ),o C 3+2,o) 



(D.6) 
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( c ]+i,ct c j> ~ c ],a £ !j+i )C r)(nj_i,_ ( j + n.^-^ + % + i _ CT + n j+2 - a - 2)] 



- n 



"3-o 



+ ( C ]+l,* C 3,° + c ],<T C 3+l,°) (4 ,-a c 3-l -o + C )-1-<t C 3 ,-<0 ( n j+l,-^ ~ ^) + 
+ (4+2,-a C i+l,-^ + C ]+l ,-^ C J+2 -a) (nj-a - ^ } 



(D.7) 



These operators have the following properties, (i) they are Hermitian and renormalized, 
22, [0) = 0, (ii) they are SO(4) invariant, [H s ,S a ] = = [H s ,rj a ], (iii) 22 S (C7 = 0) = F s , (iv) 
their action on the one-particle states \k,a) = cJ.(A;)|0) is obtained as 



Hi\k,a) 



-2 cos k — — 
2 



\k,a), 



H 2 \k,a) = [2sin2/c + 2£/sinfc] \k,a), 



H 3 \k,a) 
H 4 \k,a) 



3U 

2 cos 3k — h 3U cos 2k + U 2 cos k 



\k,a) 



3U 2 U 3 
-2 sin 4A; + U{4 sin A; — 4 sin 3A;) — sin 2A; -\ — — sin A; 



(D.8) 
(D.9) 
(D.10) 

cr>. (D.ll) 



The latter formulae lead us to the conjecture ( |5.24| ) in section 5.3. We make a remark here 
that there is an arbitrariness in the definition of higher conserved quantities by adding linear 
combinations of lower ones. We do not know what is the most natural choice for higher 
conserved quantities. 
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